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We determine the response of a self-similar isothermal stellar system to small adi- 
abatic gravitational perturbations. For odd spherical harmonics, the response is 
identical to the response of the analogous isothermal fluid system. For even spher- 
ical harmonics, the response can be regarded as an infinite series of wavetrains 
in logr, implying alternating compression and rarefaction in equal logarithmic 
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t""-. , fields from external sources are not strongly amplified by an intervening isother- 

mal stellar system, except at radii < 10~ 3 ' 5 times the satellite radius; at some 
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radii the stellar system can even screen the external tidal field in a manner anal- 
ogous to Debye screening. As Weinberg has pointed out, individual resonances 
in a stellar system can strongly amplify external tidal fields over a limited radial 
range, but we cannot address this possibility because we examine only adiabatic 
perturbations. We also discuss the application of our method to the halo response 



caused by the slow growth of an embedded thin disk. 



Subject headings: stellar dynamics - galaxies: individual (Milky Way) - galaxies: 
haloes - galaxies: kinematics and dynamics 



1. Introduction 

In models of galaxy formation based on hierarchical clustering, the gravitational field be- 
comes less and less smooth at larger and larger distances from the galaxy center. The halos 
of isolated galaxies contain satellite and companion galaxies and merging dark matter sub- 
structure, and galaxies in groups and clusters are subject to slowly varying tidal fields from 
distant group members as well as rapidly changing forces from close encounters. Despite 
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this noisy environment, the visible inner regions of most galaxies are relatively smooth and 
regular. 

One expects that gravitational noise from the outer halo and beyond will not severely 
disturb the inner galaxy because the dominant (quadrupole) tidal potential at radius r from 
a source of mass m p at radius r p ^> r varies as Gm p r 2 /r|; in other words the ratio of the 
tidal force to the force from the body of the galaxy itself is of order (m p /m)(r /r p ) 3 where 
m is the mass of the galaxy within r. The cubic dependence of the force ratio on the radius 
ratio implies that tidal forces from distant satellites are generally unable to strongly perturb 
the inner galaxy (r <C r p ). 

However, this assessment neglects an important effect: in a realistic galaxy model the 
tidal forces do not propagate through a vacuum — rather, they propagate through the density 
field of the dark-matter halo. The halo provides an intervening medium that modifies the 
external tidal field in a manner analogous to the polarization of a dielectric medium in 
electrostatics (Jackson 1975). Another closely related analogue is the dimensionless Love 
number in geophysics (e.g. Jeffreys 1970), which measures the ratio of the direct tidal 
potential from the Moon to the augmented tidal potential that includes the gravitational 
potential arising from the deformation of the Earth in response to the lunar tide. 

The influence of the halo response on tidal fields was first pointed out by Lynden-Bell 
(1985), who computed the Love number for a spherical fluid halo with power-law density 
distribution, po(r) oc r~ a (see also Nelson & Tremaine 1995); we shall find below, however, 
that the response of fluid and stellar systems with the same density distribution can be 
quite different. Weinberg (1995, 1997) found that the response of the Galactic halo strongly 
enhances the direct tidal field from the Magellanic Clouds and suggested that the resulting 
disk distortion could account for the location, position angle, and sign of the HI warp in the 
outer Galaxy. In many respects Weinberg's calculations are much more sophisticated than 
ours, as we shall focus on a simplified case that provides insight into the phenomenon rather 
than accurate numbers for a realistic system. 

In the present paper, we examine the response of a spherical, scale-free isothermal stellar 
system to an adiabatically applied external gravitational perturbation. This problem only 
approximates reality in that it neglects the time- dependence that accompanies most tidal 
fields (e.g. those due to orbiting satellites). The reward is an analytically tractable problem 
that permits us to understand the physics of the halo response — and the linear response 
of stellar systems in general — in a more detailed manner than has hitherto been possible. 
Moreover, we are most interested in the response of the galaxy at small radii, where the 
dynamical time is short and the approximation of a static tide is not unreasonable. 
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Adiabatically applied gravitational perturbations are also relevant to other astrophysical 
phenomena, including the slow growth of a central black hole in a galaxy (e.g. Young 1980, 
Quinlan et al. 1995) and the slow growth of a galaxy disk in a halo or spheroid (Binney & 
May 1986; Dubinski 1994). Thus we shall also apply the tools we have developed to these 
problems. 

The plan of the paper is as follows. In §|2|, we develop the solution to the linearized 
Boltzmann-Poisson equation that describes the response of a stellar system to an adia- 
batically applied perturbation. The relationship between adiabaticity and reversibility is 
discussed and a comparison with analogous fluid systems is also made. In §|3|, we derive 
the adiabatic Green's function for the response of the isothermal sphere. This machinery is 
applied in §[| and the results are discussed in §|5|. 



2. Linearized Adiabatic Response 

We examine the response of a stationary stellar system with unperturbed distribution func- 
tion (hereafter DF) F(r,v) and gravitational potential $o( r ) to a weak external potential 
$ e (r, t). The evolution of the perturbed DF /(r, v, t) is described by the linearized collision- 
less Boltzmann equation (e.g. Kalnajs 1971; Weinberg 1989), 

df df df „^ df „^ dF 

— = — + v • — — V$ " — = V$i ■ ; (1) 

dt dt dr dv dv 

the first equality defines the convective or Lagrangian derivative in phase space. The per- 
turbing potential $i(r, t) = $ e (r, t) + <& s (r, t), where $ s is the response potential, given by 
Poisson's equation 

V 2 $ s = 4ttG J fdv. (2) 

We shall restrict ourselves to the case where the DF of the stationary stellar system 
depends only on energy, F(r, v) = F(E), where E = |v 2 + $o( r ) is the energy per unit 
mass. Then equation (|1|) becomes 

d f dF 17 

- = v.V 9i — = Fe \ — -—.) ) ( 3 ) 

where Fe = dF/dE. Note that — v • V$i is the power per unit mass delivered by the 
perturbing potential. 



In this paper we shall focus on the case where the external perturbation grows slowly 
from zero in the distant past. We assume that the potential < &o( r ) admits action- angle 
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variables (I, w), so that the unperturbed stellar orbits have the form I =constant, w = 
u;(I)t+constant. We then expand the perturbing potential in a Fourier series of the form 
$1 = Ylic ^k(I) t) exp(ik • w), where k is an integer triplet. We assume that 3>k(I, £) = 
$k(I) exp(et) where e is small and positive. Then we can integrate (ffi) to obtain 



/ = F £ 5> k (M) 



k- 



k- 



w — le 



exp(zk • w). 



(4) 



In the limit e — > the square bracket approaches unity, unless k ■ w = 0, in which case it is 
zero. The condition k ■ w = can be satisfied in three different ways: 



1. At isolated or local resonances in phase space, where k • w(I) = for a particular set 
of actions I (accidental degeneracies). We shall ignore such resonances since they are 
unlikely to dominate the overall response of the stellar system. 

2. When k = 0. 

3. When some symmetry of the stellar system dictates that there are non-zero values of 
k such that k • a>(I) = for all values of the actions I (global resonance). For example, 
with a suitable choice of actions, in spherical potentials u 3 = (the orbital plane does 
not precess) so that k • u> = whenever k\ = k 2 = 0; and for Kepler potentials u\ = u>2 
so that k ■ u) = whenever k\ = — k 2 . 



For a given stellar system, we shall denote by N the set of integer triples k such that 
k ■ u>(T) = for all values of the actions I (i.e. all k satisfying conditions 2 or 3 above). Then 
in the adiabatic limit e — > equation (^) becomes (e.g. Lynden-Bell 1969) 

f = F E Y, $ k (I,t)exp(ik- w) = F E [& 1 - ($!>], (5) 

where = J2k&N $k(I, t) exp(ik • w). can be regarded as a time average of the 

perturbing potential, over times long compared to the orbital period but short compared 
to e _1 . This is not the same as an average over orbital phase, because the time average of 
exp(ik ■ w) 7^ if k £ iV but the phase average of exp(ik • w) ^ only if k = 0. 

Note that the solution (^) is only one of many possible static solutions to the linearized 
collisionless Boltzmann equation if the set iV is not empty. If we write the perturbed DF in 
action-angle variables, / = Z)k/k(I) exp(zk • w), then the linearized collisionless Boltzmann 
equation reads 

£exp(2k-w)(k-u,)(/ k -F £ $ k ) = 0, (6) 

k 
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which has the general solution 

/ = F E $ 1 + Yl exp(ik • w) (7) 

where Pk(I) is an arbitrary function; equation @ corresponds to the particular case pk = 
-F E $ k . 

The absence of the term k = in equation (|5|) has a simple physical interpretation in 
terms of entropy. The entropy of the unperturbed stellar system is 

Si = J drdvF In F = (2vr) 3 J dIF In F, (8) 

because F(E) is a function of the actions only and drdv = dldw. The final entropy is 

S f = J dIdw{F + /) ln(F + f) = Si + J dldw(l + InF)/ + 0(/ 2 ). (9) 

Using equation @, we find the change in entropy 

AS = S f - Si = f dItZw(l + \nF)F E ^ $ k exp(zk ■ w), (10) 

which vanishes upon integrating over w, because OGiV. 

We can also compare the solution (|5|) to the response of a barotropic fluid system to 
a weak, slowly growing external potential. In an equilibrium barotropic fluid, the pressure 
P, density p and potential $ form a one-parameter family under an isentropic change; thus 
P = P(p) = P(&)- The equation of hydrostatic equilibrium states that dP = —pdQ so 

dP , s 

For a small change in potential $i, we can determine the perturbed density by a Taylor 
expansion of this equation. Defining p = p + p 1; we find that 

For comparison to stellar systems, take equation (^|) with pk = 0, and integrate over velocity: 

pi = | fdv = $i | i^civ. (13) 
Since the unperturbed density po — j F(\v 2 + $o)^ v > we have (dp/d§) = J F E dv, so that 



o 
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The result (|14]) is equivalent to flI2|) . In other words the response of a barotropic fluid system 
to a slowly growing potential is the same as the static solution 

fi = Fe*! (15) 

of the stellar system with the same density distribution, but is not the same as the response 
of this stellar system to the same slowly growing potential. 

In the following sections we shall examine the adiabatic response of both stellar (eq. ^) 
and fluid (eq. [15]) systems; the fluid is less realistic but captures much of the relevant physics 
with less algebraic complexity. 



3. Adiabatic response of an isothermal sphere 

We now investigate the solutions to equations (0) and (|) in the case where the stationary 
stellar system is a singular isothermal sphere (Chandrasekhar 1939, Binney & Tremaine 
1987), in which the DF F(E), the density po(r) = J Fdv, and the potential $o( r ) have the 
form 

F(E) = Kexp(-f3E), p (r) = ^f, <£> (r) = |log(r/r a ), (16) 

r z p 

where 2%Gp a r 2 a p = 1, K = (2n)-^ 2 p a f3* /2 . Equations (@) and (§ can be written as 

V 2 <T + ^($ s + $ e ) = -AnG J dvF E (<& s + $ e ) r>v , (17) 

where (-) r ,v denotes an average — in the sense of equation ([5]) — over the orbit passing through 
(r, v). This constitutes an integro-differential equation for the response potential $ s . In the 
case of a fluid the right side of equation (|T7|) would be zero. 

Write 

®(r) = J2®t™(r)Y em (n). (18) 

tm 

Then after multiplying fll7|) by Y^ m (fl) and integrating over dfl we get 

-AnG J dQY/jn) J dvF E (T, (*J m , + n, m ,) Y e J) • (19) 

\i'm' ' r v 

We now multiply ([OJ) by r a+1 and integrate over r. We assume that $| m ~ r~ a ~ as r — ► 
and ~ r ~ a + as r _ > 00; where a+ < £ + 1, a_ > — £ since otherwise the density distribution 



-7- 



that gives rise to $| m is unphysical. We restrict a to the strip of the complex plane defined 
by 

-i < a_ < Re(a) < a+ < £ + 1, (20) 
so that boundary terms arising from integration by parts vanish. 

The Mellin transform of a function y(r) is written y(a), where 

roo | ra+ioo 

y{a) = / r Q_1 y(r)dr, y(r) = — / y(a)r~ a da. (21) 

JO 27T? Ja-ioo 

Thus 

a(a - 1)$L(«) + [2 - + l)fe(a) + 2$? m (a) 

= -4ttG / drdvr*~%* m (n)F E (]T(<t>|, m , + ^w)w\ . (22) 

We now convert to action-angle variables (I, w). Our conventions follow Tremaine & 
Weinberg (1984): I 2 is the total angular momentum, J 3 is the z-component of the angular 
momentum, and I\ is the radial action, 

1 t Ta 

h = - v r dr, (23) 



7T Jr. 

where v r = [2E — 2$ ( r ) — ^l/ 7 " 2 ] 1 ^ 2 * s the radial velocity and r p and r a are the pericenter and 
apocenter distances, at which v r — 0. The only angle with a direct geometrical interpretation 
is u>3, which is the azimuth at which the orbit crosses the equatorial plane upward (the 
ascending node). Motion in the unperturbed Hamiltonian H (l, w) = |t> 2 + $o( r ) is given 
by I =constant, w = u>t + w where u> = dH /dl. Note that u; 3 = 0, since the orbital plane 
does not precess in a spherical potential, so that k ■ u = if ki = k 2 = 0. Note also that 
drdv = dldw. 

Expanding the potentials in action-angle variables, we have 
mUr') + n, m ,{r')\YE<m>{n') = £ fym'Vtwfflw'foQ!) exp (i £ l' k w' k ) , 

I'^l'z \ k=l ) 

r^-'YUn) = Y, 6 h m Vi h iMUl(l)exp(tj2hw k ). (24) 

hl 2 l 3 V k=l ) 

Here cos/3 = /3//2, ^ 2 i 3 (/3) is defined in terms of rotation matrices by Tremaine & Weinberg 
(1984), and 

= ^|^^iexp(-^ lWl )[^ m (r) + $, e m (r)]expH/ 2X (I,^ 1 )], 
U\\ 2 {I) = — f T ^ 1 exp(-?/ lWl )r Q *- 1 expH/ 2 x(I,u'i)]. (25) 
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In this equation 

rr(l,wi) d r 

X(l,w 1 )= -(u 2 -I 2 /r 2 ). (26) 

Jr p (I) V r 

We now find the average (-) r ,v of the first of equations (0). This is obtained by replacing 
I' by I, f3' by /?, w' k by Wu and restricting the summation to 1' 2 , 1' 3 ) G N. Now substitute this 
result into equation (JZ2"|), replacing drdv by dldw and r a_1 Y£ m (f2) by the complex conjugate 
of the second of equations (p4f). Only terms with ^ = l[, l 2 = 1' 2 , and Z 3 = 1' 3 = m = m' 
survive the integration over w. Moreover, in the singular isothermal sphere the set of resonant 
triplets N is given by h = l 2 = 0, since the only global resonance is due to spherical symmetry. 
Thus, the right side of (p2|) becomes 

-4ttG / dvdvr^Y^FE ( £(SJ m , + n>m>)Yt>J) 

J \i'm' I r ,v 

-2 5 tt 4 G£ I dhhdhd cos pv; 0m (P)Ul*(I)F E V il0m (P)WX m (l)- (27) 
t> J 



We next use the orthogonality relation (Edmonds 1960) 



d cos (3V; i2m (f3)V m2m (f3) 



2£+l 



in particular 



a 



i 



CO 



2tt 2 
0. 



1M 2 

2 ' 



r(i£ + i) 



even. 



odd. 



(28) 



(29) 



Thus equation (p2|) becomes 

q(q - i)$L(«) + [2 - + i)]$L(«) + 2$L(«) 

= -2 5 n 4 C eo G J dhI 2 dI 2 F E U°*(I)W° 0m (I) 

= -8n 2 C m G J dIJ 2 dI 2 F E J dw 1 r a -\l,w 1 ) J dw[ ($| m + $, e J [r(I,w[)} 

RdRdEIl (E) F E f ^ y ^ ($L + $|J [r(I> ^ 



(30) 



in the last line we have changed the integration variables from I\,I 2 to the energy E and 
the dimensionless angular momentum R = I 2 /I C (E), where I C {E) is the angular momentum 
of a circular orbit of energy E. 
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Since the potential is scale- free the radius can be written in the form r(I, Wi) = r c (E) x 
x(R, Wi), where x is a dimensionless function that can be determined by numerical integration 
of the orbits. Changing the integration variable from E to r' = r c (E)x(w[, R), the right side 
of (BUI) becomes 



-8n 2 C eo G J RdR J dw l x a ~ 1 {R,w l ) 



x / dw[x- a (R, w[) I dr'r' - 1 [<&| m (r') + $f m (r')] 



dEI 2 (E)F E 



dr c u>i(E, R) 



For the singular isothermal sphere 
dEI 2 c (E)F E 1 



dr^E^R) 2 1 /H 5 / 2 eGg{R)' 
is a function only of R. Thus 



where g(R) = Ux(E, R)r c (E)f3 



1/2 



(31) 



(32) 



a(a - l)$J m (a) + [2 - £(£ + l)]$| TO (a) + 2d>| m (a 
2 5 / 2 Q fffc/fi 



7rV2 e 7 g(R) 



dw 1 x a -\R,w 1 ) dw' x x- a {R,w' x ) $| m (a) + $ £ e m («) 



where 



Note that 



iffcO 



2 5 / 2 /■ itai? 



H*(a) = H(a*), H(l-a) = H(a). 
It can also be shown after some algebra that 

if (0) = if (1) = 4vr; if' (1) = -if' (0) = tt. 



(33) 

(34) 
(35) 

(36) 



Equation fl3"3] ) can be rewritten as 

i(«) = ^^(«), (37) 

where 

D e (a) = a(a - 1) + [2 - l{l + 1)] - C eo H(a). (38) 
Both if (a) and D#(a) are even functions of a — |. 

Thus the linear response of the singular isothermal sphere to adiabatic perturbations 
for odd I is analytic and for even i requires only the numerical evaluation of a single function 
H(a). 
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3.1. Zero-frequency normal modes 

Roots of the dispersion relation De(a) = may be regarded as zero-frequency normal modes 
of the singular isothermal sphere. The contribution to $| m (r) from each root of Dg(a) is 
proportional to r~ a , which follows from the inverse Mellin transform. In general these 
solutions should not be interpreted as physical normal modes, since a power-law potential 
perturbation is always nonlinear at large or small radii, except in the special case Re (a) = 0. 
In the fluid case, this limitation is highlighted by Lichtenstein's theorem, which states that 
static, isolated, equilibrium fluid systems must be spherical, i.e. there are no finite normal 
modes with zero frequency (e.g. Lindblom 1992). 

Two trivial normal modes arise from gauge transformations: the mode a = 0, £ = 
corresponds to a shift in the zero-point of the potential, and a = 1, £ = 1 corresponds to a 
uniform translation of the unperturbed stellar system. 

The symmetries (|35| ) imply that if a is a root of the dispersion relation, then so is a* 
and 1 — a. Thus if x = a— | the roots come in pairs \-\-x, \ — x when x is real or imaginary, 
and in quartets | + x, | — x, | + x*, | — x* if x is complex. 

Finding the roots of Di(a) is easy when £ is odd: then Ci = so that the roots are 

a = \ ± [£(£ + 1) - |] 1/2 , £odd. (39) 

To find the roots when £ is even, it is convenient to rewrite the function H(a) explicitly as 
a complex function of a complex variable. Defining a = a + bi where a and b are real, we 
find the real and imaginary parts of H(a) to be 

Re\H(a)] = / ^f^- [ dw 1 [ dw' l cos[bln(x/x')]x a ^ 1 x 

7r i / i e J g{R) J J 

o5/2 „ pj D /. /. 

Im\H(a)} = / ^ / dwx / dw[ sm\b\ii{x/x')}x a - 1 x'- a . (40) 

7r i / i e J g{R) J J 

The roots of Di(a) occur at the simultaneous zeros of its real and imaginary parts. 
Figure 1 shows the zero curves in the complex plane of the real and imaginary parts of 
Dt(ot), for several values of £. The roots with the smallest values of \a\ = | Re(a)\, which we 
denote a = ±i6 and ±ibo, have the greatest physical significance, since these generally 
determine the asymptotic behaviour of the response of the system at small and large radii. 
These roots are given in Table |l|; note that they satisfy the constraint — £ < a_, a + < £ + 1 
given by equation (|20D . Also note that + clq = 1. 



a— 1™'— o 



The normal mode corresponding to the root Oq = 1 for £ = 1 is simply a uniform 
displacement of the unperturbed system. 




Fig. 1. — The zero curves of the real (dotted) and imaginary (solid) parts of D^a) (eq. [38]) 
for indicated values of £. The roots (solid points) lie at the intersections of the two curves. 
Those with largest real parts in the left half-plane will dominate the response at small radii. 
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For comparison, the dispersion relation in the fluid system is obtained by setting H(a) = 
in equation fl38|). The corresponding roots are given by ( |39| ) for all £, both even and odd. 
Thus for i = the fluid system has a = |(1 ± iyf); this is a familiar result since the 
difference in potential between the singular and non-singular isothermal spheres oscillates 
with this spatial frequency at large radii (Chandrasekhar 1939). For £ = 2, equation ([39]) 
yields a = \{l ± y/Tf). 



3.2. Green's function 



We now discuss the adiabatic response of the singular isothermal sphere to a point mass m p 
at position r p , which corresponds to the external potential 

nr)=Y,®Ur)YUn), where $| m (r) = -— -f Yl m {n p )-^. (41) 

Here fl p denotes the angular coordinates of r p , r < = min(|r|, |r p |) and r> = max(|r|, |r p |). 
The Mellin transform of the point-mass potential is therefore given by 

AnGrrip x ( 1 1 



= -^fYLW {— e ~ ^J^j) • -I < < 1+ I- (42) 

We now employ equation (p7|) to find the response 



Table 1: Asymptotically dominant roots 





a 


aS 


bo 





0.0000 


1.0000 


0.0000 


1 


0.0000 


1.0000 


0.0000 


2 


-1.5623 


2.5623 


0.4804 


3 


-2.7016 


3.7016 


0.0000 


4 


-2.0460 


3.0460 


0.5480 


5 


-4.8151 


5.8151 


0.0000 


6 


-2.2771 


3.2771 


0.5457 


7 


-6.8655 


7.8655 


0.0000 
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The inverse Mellin transform is 

2m Ja-ioo 



$| m (a)r a da 



Gm p 2i Y * (ft ) 



r 



!> 



a+ioo 



da 



a + 



a — 



C m H(a) 



D P ( a) 



(44) 



The constant a must be chosen so that — £ < a < £ + 1 and a_ < a < a + (eq. |20[). If we 
choose a = \ and use the fact that H(a) and Dp(a) are even functions of a — ^ then it is 



easy to show that 




where 



^ lm {u) = #/ m (l/ii); 



in other words the response for r < r p is determined by the response for r > r p . 



(45) 



To evaluate the integral ([44]) in closed form, we analytically continue the integrand over 
the whole complex plane and then close the integration contour, in the left half-plane for 
r < r p and in the right half-plane for r > r p . The contributions to the integral come from 
the poles at a = —£, a = £ + 1, and the roots of Dp(a), which we denote by a = | ± Xk, 
where Re(xk) > 0. The asymptotic behavior at small radii, ~ r~~°~ , is determined by 
the pole in the negative half-plane with the smallest | Re(a)\ (either from Table |l]or —£), 
while the asymptotic behavior at large radii, ~ r~ a+ , is determined by the pole in the 
positive half-plane with the smallest | Re(a)\ (either ciq from Table p] or £ + 1). Thus 



AirGm 
W+T)r p 

AnGm 




r < r 



p- 



(2£ + l)r 



' p 

r , 



J2£+]_ 



r > r 



p- 



(46) 



these two expressions satisfy the symmetry relation fl45|) . The sum is over all roots with 
Re(a — |) > 0; recall that if Im(a) ^ then a* is also a root. 

The first term on the right side of each equality defines a response that precisely cancels 
the external perturbing potential $£ m (eq. [II]). The total potential perturbation §\ m = 
$| m + therefore is entirely determined by the spectrum of normal modes. When i is 
even the roots x^ are generally complex (Table P so ^\ m {f) varies sinusoidally in logr; for 
odd i the total potential has no phase variation because the roots are real. 

Since D[(a) is an odd function of a — |, the response potential $| m (r) is continuous 
through the shell r = r p . The total response $ s (r) = J2e, m ^ ) lm( r )^m(^) is also continuous. 
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In addition, d§\ m jdr is continuous through the shell r = r p ; this can be shown using the 
identity 



V - 



u m+^) (47) 

which in turn can be derived by considering the integral J c da(a — \)/Dp(a) where C is the 
circle lal — > oo. 



In the case of a fluid system, the analogous expressions are 



AnGm 



p v-* 



{2£ + l)r/ /m( " p) 



r 



2£ + 1 / r \ 
1) - 71V2 l^j 



[£(£+l)-7/4] 1 / 2 -l/2' 



47rGm 



(2£+l)r p 



e+i 



[4£(£+l)-7] 

21 + 1 
[A£{£ + 1) - 7] 1 / 2 



[^(^+l)-7/4] 1 / 2 +l/2' 



, r <r p 



, r > r p . 



(4f 



A similar result was derived already by Lynden-Bell (1985), with minor differences — in par- 
ticular he assumed that the fluid halo stopped at r p , which changes $£ m by a multiplicative 
factor. 

Some care is required to interpret these expressions when £ = or £ = 1. For £ = 0, the 
response potential interior to the point mass contains constant terms (oc r e and oc r Xk ~ 1 ^ 2 for 
x k — §)i these terms exert no force and can be eliminated by re-defining the zero-point of the 
potential. When £ = 1, the concept of the linear adiabatic response to a small perturbing 
potential is ill-defined, since a small but steady perturbing potential can lead to a large 
shift in the center of mass of the stellar system. This indeterminacy can be reflected in 
equations ( P>D and (^) by adding an arbitrary amount of the normal mode oc r _1 that 
corresponds to a uniform translation of the unperturbed stellar system. 

In practice, it is easiest to determine the response potential through direct numerical 



calculation of the inverse Mellin transform ([44]) , rather than by summing the residues at all 
the poles. The simplest choice of contour is along the symmetry axis, a = |, for which the 
imaginary part of H(a) vanishes. The real part of equation 



then becomes 



rr/l ,N 2 5 / 2 [RdR\ 



dwi cos(Mnx)x 



-1/2 



dwi sin(61nx)x" 



-1/2 



(49) 



which has a narrow peak about 6 = and decays approximately as b 1 for b > 50. 

It is also helpful to rewrite equation ( |44] ) as a Fourier transform which can be evaluated 
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using FFTs (e.g. Acton 1990): 

= -g^yUV) f *«» V-Hr/r,)] u % I it) ' <50) 

The integration parameter b represents the logarithmic wavenumber of the response. In the 
fluid case, 

AGtyi f°° 2 

*'^ r) = -J^^W* Jo ^ CQS[Mn(r/rp)] [( i + £) 2 + &2 ][(i + £) 2 + &2 -2] - (51) 



3.3. Density response 

The density response is determined directly from the Green's function (f46|) through Poisson's 
equation. Write the density as a multipole expansion: 

p s (r) = J2ptm(r)Y im (n); (52) 
then substitute into Poisson's equation so that 

^| m (r) = ~/ d ^ - ^±il*J m . (53) 

There is no surface- density layer at r = r p , because the response potential and its gradient 
are continuous at r = r p (see the discussion following equation |56|). 



Substituting equation ([46]) for $| m gives 



l\2l 



r l k\ r pJ D A- 2 - x k) 

- m *Y* f o )V ^V fc+5/2 [xi " (£ + ^ )2] r>r 



(54) 



At small radii, the fractional density perturbation is p\ m j po ~ r a o , while at large radii the 
fractional perturbation varies as r~ a o~ (cf. Table [I]). 

For a fluid system, 

m 2Y* (fl ) ( r \ W+V-W 1 ' 3 - 5 / 2 
PlmK ) ~ i% [U{1 + 1) - 7]V2 \ rp } ' S r P 



2y; m (^ p ) /r r ,\ [^(^+i)-7/4] i / 2 +5/2 



r3 [4£(£+l)-7] 1 / 2 Vr 



r > r p . 



(55) 
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3.4. Amplification factor 



It is useful to express these results in terms of the overall amplification of the point-mass 
perturbation by the response. Defining the amplification Xi( r ) — (®e m + ^D/^L) we fi n d 
that 



The amplification also obeys a symmetry relation similar to (fl5|), xt( r ) — Vv( r / r p) where 
tp e (u) = ipe(^/u)- 

As we discussed after equation ([18|) these formulae are not meaningful for I = when 
r < r p , or for 1=1. 

Asymptotically Xe( r ) ~ r^-?- 1 / 2 — r -% ~ l asr o, and Xi( r ) ~ r ~ x o+t+i/2 _ r -% +e+i 
as r — > oo. Values of ao are given in Table |1] and we see that for all i > the amplification 
diverges as r — > or r — > oo. However, for even £ this asymptotic behavior does not 
necessarily appear until very large values of | log r | are reached. For odd £, there is no 
distinction between the exact and asymptotic behavior since the response is determined by 
a single root. 



Consider the response of the halo to a point-mass satellite on a circular orbit of radius r p . 
Our assumption that the satellite perturbation is adiabatic is not accurate at radii r ~ r p , 
where the characteristic orbital period in the halo is comparable to the satellite's orbital 
period. Nevertheless, the adiabatic approximation is plausible at radii r <C r p , and should 
approximately describe how the tidal field of the satellite is propagated to small radii. The 
monopole [I = 0) response of the halo is not so interesting, since it is difficult to distinguish 
observationally from the potential of the unperturbed stellar system. As we have discussed, 
our formalism is not powerful enough to determine the dipole (£ = 1) response. Thus we 
shall focus on the quadrupole (£ = 2) response. 

Figure 2 shows the amplification X2(f ) of the quadrupole tidal field from a point mass at 
radius r p . The amplification has a weak maximum of about 2 at log(r/r p ) m —1 and becomes 




(56) 



4. Applications 



4.1. Tidal amplification 
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negative for log r/r p & —2. This oscillation arises because the asymptotically dominant root 
in Table |I] is complex. The amplification remains negative down to log r/r p pa 5 — far too 
small to be of interest — where it begins a positive rise. The dashed line shows the much 
stronger quadrupole amplification for the analogous fluid system. 

These results indicate that the response of an isothermal halo does not greatly enhance 
the quadrupole tidal field from an external point mass: the response is much smaller than 
for an isothermal fluid and over a large range in log(r/r p ) the halo response actually screens 
or suppresses the external tidal field. 



4.2. Response to disk growth 

We can also calculate the response of the halo to the adiabatic growth of an embedded 
disk (Binney & May 1986; Dubinski 1994). We compute the response using scale-free, 
axisymmetric, razor-thin disk models, which have density distributions 

/ x x <$(cos#) /K _s 
pr =£r- 57 

r 

We shall use the relation 

5(cos9) = 27r5^io(|M)*«,o(0,O); (58) 

i 

note that only even harmonics contribute to the sum. 

Using equation (^), we calculate the total potential arising from the embedded disk: 

i 

1 r°° , , ( r\ Xk ~2 



/V(rT +2 E(r')}. (59) 



D' 2i (l + x k ) Jo 



For example, in the case of a Mestel disk, S(r) = S a r a /r, the monopole potential is 

, , . 4vr 3 / 2 GS a r ai . . 

^oo(r)= jy {1) — Inr; (60) 

here we have neglected constant contributions to the potential (including a divergent one) 
which change the zero-point but have no physical effects. For i > 0, we have 

$^ (r) = 167T 2 GS a r a F 2 * j0 (±7r, o) £ „ (1 -gT - (61) 
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Fig. 2. — Amplification X2{ r ) of the quadrupole tidal field from a point mass at r p . The 
solid line shows the exact result from numerical evaluation of equation fl56|). The dotted line 
that follows the solid line shows the amplification due to the asymptotically dominant root. 
The dotted line that oscillates weakly about zero shows the contribution from the next most 
important root. The dashed line shows the much stronger amplification for the fluid case. 
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These terms give rise to tangential forces but no radial forces. 

We can evaluate sums over roots by analogy with equation (fT7|), using the integral 
f c da(a — 2) / Di(at)(at + 1 — 7) (a — 2 + 7). This leads to the identity 

?^(|-^)[(|-7) 2 -4] = 2^(2- 7 )' (62) 
Using this result for 7 = 1 and the identity (|36|), equation ( |61] ) simplifies to 

f , _ 87r 2 GS a r a ^ (tt/2, 0) 
$2 ^ o(rj " 2-2*(2*+1)-4ttCV ( j 

which, remarkably, is completely analytic. 

We can easily extend these results to disks with surface density E(r) = S a (r a /r) 7 , 
1 < 7 < 2. The total potential is 

W) = ^ (2 _ 7) (7 J ^,0 0) . (64) 

We can also express this result in terms of the amplification \2i — ^21 0/^21 0' w here <& d is 
the direct potential from the disk. We have 

= (I-7) (2 -7) -21(21+1) 

X2e (i_ 7 )(2- 7 )-2^+l) + 2-C^(2-7)- 1 J 

This expression is valid for all values of 2£ and 7 G [1,2] except for 2£ = 0, 7 = 1. Note that 
for large £, \2t 1; the halo does not amplify the response from the disk on small scales. 
Table shows values of X2io for a range of 2£ and 7. 



Table 2: Disk amplification 



2£/ 7 


1 


5 
4 


3 
2 


7 
4 





2.00 


1.93 


1.90 


1.93 


2 


1.39 


1.33 


1.33 


1.33 


4 


1.10 


1.10 


1.10 


1.10 


6 


1.05 


1.05 


1.05 


1.05 


8 


1.03 


1.03 


1.03 


1.03 



Using Poisson's equation, we derive the response density of the halo in the presence of 
the disk. The response consists of a density 

P S = ^? " ^o&r, ( 66 ) 
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This can be expressed as an induced enhancement S in the local halo density 

S = ^= 2 ffX - 1)1^0(5^ °) W> 0), (67) 

P0 Pa' j> 

where p a is defined in equation flip]). Figure 3 shows contours of the fractional density 
enhancement 5 for various values of 7. 



4.3. Response to black hole growth 

We may also consider the response of the halo to the adiabatic growth of a central dark 
object, such as a massive black hole. In this case the nonlinear problem has already been 
investigated and solved numerically by several authors (Peebles 1972; Young 1980; Quinlan 
et al. 1995), so the results of our linear calculation are only of academic interest. 

We take the limit of the second of equations fl46|) in which r p —>■ 0, and consider the 
monopole case £ = m = 0. Only the term corresponding to the dominant root | + — 1 
survives since $| m oc r p k ~ 1 / 2 and all other terms have Re(xk) > \- The total potential is 
therefore 

= ($00 + $00) ^oo(tt) = ~ L W)Y l , (68) 

where D' (i) = ~ so that the linear response amplifies the direct black- hole potential by a 
factor of 2. 

The interpretation of this result is as follows. At radii r < = /3Gm p the black hole 
induces a (nonlinear) density cusp in the stellar system. At radii r > rj, we might expect 
a linear density response, but, as it turns out, the linear terms vanish so the response is 
second-order in the black-hole mass (this is straightforward to show in the case of a system 
composed of stars on circular orbits). Thus the total potential perturbation at large radii 
consists of the Keplerian potential of the black hole augmented by the mass in the density 
cusp. The particular value of the augmented mass given above, m p /D' (l), is not necessarily 
accurate since it is derived from a linear calculation while the cusp is nonlinear. 



5. Discussion 

We have discussed the response of a self-similar stellar system to weak adiabatic gravitational 
perturbations. Our problem is idealized because real stellar systems are only approximately 
self-similar, and because we neglect the time-dependence that accompanies most tidal fields, 
such as those from orbiting satellites. Nevertheless, this calculation provides one of the very 
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Fig. 3. — Contours of the fractional halo response density 5 (eq. ^) induced by razor-thin 
disks with surface density S a (r a /r) 7 . Contours show levels 5 = 0.2,0.3, . . . ,0.8 in units of 
^a/faPa, starting from the topmost contour in each panel and increasing toward the disk. 
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few examples where the linear response of an inhomogeneous stellar system can be explicitly 
computed in terms of quadratures, and hence offers analytic insight into the nature of the 
response. Our results can also be used to test numerical codes used in linear response 
calculations. 

Our assumptions that the stellar system is self-similar and that the perturbation is 
adiabatic can be lifted in more realistic calculations using matrix methods, such as those 
of Weinberg (1995, 1997). Weinberg shows that in some cases individual resonances in the 
stellar system can lead to strong enhancements in the external tidal field over a limited range 
of radii. 

Using the methods derived above, we have investigated the propagation of the tidal 
disturbance from a static satellite into the inner galaxy; in effect, we have computed the 
static Love number for an isothermal stellar system. For odd spherical harmonics (other 
than i — 1, for which the adiabatic approximation is not self-consistent), the response 
of the stellar system is identical to the response of the barotropic fluid system with the 
same density profile. For even t > 0, we find that the zero-frequency normal modes of the 
isothermal sphere, which have the spatial dependence exp(— a logr)Y^ m (f2), have complex 
spatial eigenvalues a and hence oscillate in logr at fixed angular position Q. Partly because of 
this oscillatory behavior, the halo response does not strongly amplify the satellite's tidal field 
in the inner galaxy: the amplification factor or Love number for the dominant quadrupole 
(£ = 2) tidal component is never more than about 2 for 0.01 < r/r p < 1, crosses zero near 
r/r p ~ 0.01 and is negative over the next decade in radius (Fig. 2). 

Our calculations have some similarity to Young's (1980) investigation of the response of 
an isothermal stellar system to the adiabatic growth of a central black hole. The principal 
differences are that (i) Young's method is fully nonlinear, whereas ours is only linear; (ii) our 
method can be applied to non-spherical perturbations, such as the slow growth of a disk. 

Finally, although we have examined only the singular isothermal stellar system with 
density po( r ) ^ r <2 , h is straightforward to extend these calculations to other self-similar 
systems with density po(r) ^ r ' k - 

This research was supported by NSERC, the Fund for Astrophysical Research (CM.) 
and an Imasco Fellowship (S.T.). 
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